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ABSTRACT 



The main goal of this research is to model a flexible missile with structural 
flexibility utilizing the Equivalent Rigid Link System (ERLS) with an enhanced 
natural mode discretization. Dynamic analysis of the flexible missile in 
2- Dimension motion is presented. 

Computer simulation is performed where the pitch angle of the missile is 
controlled with a rigid-body controller. The effects of increasing payloads and 
speed to the system performance are discussed. 



THESIS DISCLAIMER 



The reader is cautioned that computer programs developed in this research 
may not have been exercised for all cases of interest. While every effort has been 
made, within the time available, to ensure that the programs are free of compu- 
tational and logic errors, they cannot be considered validated. Any application 
of these programs without additional verification is at the risk of the user. 
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I. INTRODUCTION 



A. BACKGROUND 

With the introduction of long slender missiles such as the Vanguard, the 
Redstone, and various ballistic missiles, the problem of the structural flexibility 
became severe [Ref. 1]. Due to the limited thrust available from rocket engines, 
these missiles had to be as light as possible.This meant a sacrifice in structural 
rigidity. 

Missile flexure causes additional aerodynamic loads which in turn cause ad- 
ditional flexure. Coupling occurs between the elastic modes and the control sys- 
tem as the control system gyros sense the flexure motion and the rigid body 
motion. It has become necessary to actively control the flexible structure and 
thereby reduce the structural loads and improve the vehicle response such as po- 
sition, velocity and acceleration. Reduced structural loads will also offer potential 
for reduced bending stress and fatigue problems. 

The objective of this thesis is to develop a dynamic model for a flexible missile 
and study the dynamic behavior of the flexible missile. Simulation is a valuable 
tool in the design of new missile systems and in the modification or evaluation 
of existing systems. A missile simulation allows the engineer to evaluate his design 
without the expense of actually building and flying the actual missile. System 
dynamics can be investigated through simulation with a substantial savings in 
time and expense [Ref. 2]. 

B. LITERATURE REVIEW 

Flexible missile modeling centers on the relationship between the large, rigid- 
body motion and the small motions due to structural flexibility. 

Jenkins [Ref. 2] expresses some techniques used in deriving the equations of 
motion of a rigid missile for a six degree-of-freedom (6-DOF) simulation. The 
rigid missiles are characterized by their larger size and low ratio of payload to 
total weight. Rigid missile dynamic equations were developed using the 
Newton-Euler Method. The moments and forces along with the mass and mo- 
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ments of inertia are assumed to be known in the body coordinate frame. The 
transformation between global - and local coordinate frame is achieved with a 
non-homogenous coordinate transformation matrix [Ref. 3: pp. 342]. The result- 
ing rigid-body equations of motion produced the understanding for the derivation 
of the dynamic equations of the flexible missile. 

The Equivalent Rigid Link System [Ref. 4] describes the large-motion 
kinematics of the system by the ERLS motion and the small motion relative to 
the ERLS. The application of the finite element techniques and Lagrangian dy- 
namics produces two sets of coupled, nonlinear, ordinary differential equations 
of motion, of which one set is for the large motions and the other set for the small 
motions. The small motion is described by the superposition of vibration modes. 
The modes of the vibration of the flexible bar was derived with the simple-beam 
theory [Ref. 5: pp.22 1]. In simple beam theory, it is assumed that the rotation of 
the element is insignificant compared to the vertical translation and the shear 
deformation is small relative to the bending deformation. This assumption is valid 
if the ratio between the length of the bar and its height is relatively large (more 
than 10). The set of large motion equations are non-linear in both the large and 
small motion variables while the set of small motion equations are linear in the 
small motion variable and non-linear in the large motion variables. A solution 
technique called the Sequential Integration Method [Ref. 6] was developed which 
allows efficient simulations of systems with inertia coupled motions having non- 
linear slow motion (large motion) with linear fast motion (small motion). The 
ERLS model presents a complete, efficient dynamic model able to describe large 
motion, small motion and their coupling. 

An ERLS model of a flexible spacecraft boom was developed and a computer 
simulation was performed [Ref. 7]. The equations of motion were solved using 
the Sequential Integration Method. A spatial finite discretization of the boom 
structure and the application of an assumed polynomial modal response were 
utilize^ in the approximate solution to the equations of motion. This work was 
the basis work for the modeling of the flexible missiles. 
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Ganon [Ref. 8] performed an experimental validation of the ERLS dynamic 
model in a vertical plane motion. The small motion was modeled using a shape 
matrix derived from superposition of natural modes. The agreement between the 
simulation results and the experimental results justify the application of the 
ERLS using a natural-mode shape matrix to model the dynamic response of 
flexible missile. 

C. OUTLINE 

In this study, the ERLS dynamic model is used to derive the system equations 
of motion for a flexible missile in 2-D motion. Dynamic response is predicted by 
solving the equations of motion using the Sequential Integration Method. The 
application of an assumed natural mode shape and the spatial finite element 
discretization of the missile provide an approximate solution. Computer simu- 
lation for the flexible missile is performed using MATLAB on the MACINTOSH 
computer. A rigid body controller is included in the simulation to control the 
pitch angle of the flexible missile. 

The ERLS dynamic model of the flexible missile is presented in Chapter Two. 
A rigid body controller for the flexible missile is described in Chapter Three. The 
computer simulation methodology and simulation results are presented in Chap- 
ter Four. The conclusions and recommendations for future work are presented 
in Chapter Five. 
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II. MODEL FORMULATION FOR A FLEXIBLE MISSILE 



A. KINEMATICS 

In our model, 2-D inertial reference frame, i.e., X and Y, is used for the global 
motion as shown in Fig.l. The body-fixed coordinate frame, i.e., x and y, is se- 
lected to describe the missile local motion. 




The following assumptions were made for the flexible missile: 

(1) Material density is constant throughout the body and the steel was chosen 
for modeling. 

(2) A uniform circular cross section is assumed. 

The geometric configuration parameters and material properties of the flexi- 
ble missile are listed as 
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Diameter = 0.12 meter 
Length = 4.0 meter 
Material Density = 7861.05 kg//n 3 
Young's Modulus = 2.0 x 10 u pascal 

The concept of the ERLS is applied to model the kinematics of flexible mis- 
sile. The main idea of the Equivalent Rigid Link System (Fig. 2) is to separate the 
kinematics of the flexible body into large and small motions. 




v(0) 


= 


Missile Base Deflection 


<D(0) 


= 


Missile Base Slope 


v(L) 


= 


Missile Tip Deflection 


O(L) 


= 


Missile Tip Slope 


R 


of 


Global Position Vector 
Missile Base 



= Absolute Position Vector 
of the base of ERLS 

= Local Position Vector 

= Deformation Vector 

= Equivalent Rigid Link 

= Theoretical Missile 
Position 



\ 



Figure 2. Equivalent Rigid Link System (ERLS) 
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The large motion of the missile will be represented by missile base position in 
X direction X 0 , missile base position in Y direction T 0 , and missile pitch angle 0. 
The small motions resulting from flexible motion are measured relative to the lo- 
cal coordinate frame x, y. v(0) and 0(0) are the nodal displacement and slope 
of the missile base respectively. The absolute (global) position of a point on the 
flexible missile is obtained from coordinate transformations. The global position 
(R) of any point position can be defined using ERLS in terms of a local position 
vector ( r ), a deformation vector (d), and a coordinate transformation matrix (W), 
i.e., 



R = W{r + d ) 



(2-1) 



The transformation matrix (W) between the large motion and small motion co- 
ordinate is 



1 0 0 



W = 



Xo 

Yo 



cos(d) 

sin(0) 



sin(0) 

cos(0) 



(2-2) 



and the local rigid body position vector is 






1 

x 

0 



(2-3) 



The deformation vector ( d ) is expressed in terms of a nodal displacement 
vector U and a shape function 0 as 

d = <t>U (2-4) 

where 



U = [v(0) <t>(0)] r 



(2-5) 
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The shape function <fi was derived utilizing a natural-mode superposition and 
a finite element concept. The Finite Element Method (FEM) was utilized to 
discretize the flexible body displacements and assigning the nodes. Displacements 
are for each point along the missile, a function of location and time, and it is 
necessary to discretize the deformations to obtain an ordinary differential 
equation. The natural mode shape function of a beam is used to represent the 
flexural motion of the flexible missile. Only the first two mode shapes are used. 
The flexible missile is modeled as a continuous Euler-Bernoulli free-free beam, 
neglecting shear deformation and rotary inertia effects. The nodal points are the 
base points of the flexible missile. Appendix A shows the mode shape function 
derivation. In this case, <f> is found as, 




0 0 
0 0 
a b 



(2-6) 



where a and b are defined as following, 
a = /^(C^ cos pyx + cosh /? 1 jc) + ( sin p^x + sinh p^x)) 

+ F 3 (C 2 ( cos p 2 x + cosh p 2 x) + ( sin p 2 x + sinh p 2 x)) (2 — 7) 

b = F 2 {C { { cos p { x + cosh p { x) + ( sin p { x + sinh p^x)) 

+ F 4 (C 2 ( cos p 2 x + cosh p 2 x) + ( sin p 2 x + sinh p 2 x)) (2 — 8) 

Lagrangian dynamics is used to acquire the equations of motion and the 
kinetic energy of the system will be needed in the development of the Lagrangian 
expression. The absolute velocity is listed as follows, 



R = W(T + d) + Wd 



(2-9) 
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B. KINETICS 

Lagrangian Dynamics is a systematic way to derive equations of motion for 
complex systems like flexible missiles. Lagrange's equation is written as, 



d , dKE v 8KE dPE 

dt dqi dq { dq- t 



(/ 1 ,2,...,az) 



(2-10) 



where 

K.E. = Kinetic Energy 
P.E. = Potential Energy 
q, = Generalized coordinates 
Q = Generalized forces 
n = Number of degrees of freedom 
The generalized coordinates are chosen to be, 

q = [,Y 0 r 0 e v(0) 0(0)] 7 (2-11) 



The kinetic energy of the system is defined as follows, 



K.E. = ^-\R T Rdm (2-12) 

By substituting Eq.(2-1) into Eq. (2-12), the kinetic energy is written as 

KE = ^\(f T W T VV7 + 2 r T W T W(f)U + 2r T W T W(f>U 



+ U T (j) T W T W(t)U + 2U T (f) T lV T W(f>U + U T 4) T W T W<t>U T )dm (2 - 13) 

The potential energy of the system includes the strain energy of the flexible 

S 

missile and the gravitational potential energy. 

P.E. =-j\U T r T CrUdx - \R T gdm (2 - 14) 
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where 

T = Spatial derivative of the shape function 
C = Rigidity Matrix 
g = Gravitational acceleration vector 

Appendix B includes the development of the equations of motions using 
Lagrange equations. 

Generalized forces will be found by virtual work principle. It is assumed that 
the only force is the thrust force which is applied to the base of the flexible missile 
as shown in Fig. 3. Other forces like aerodynamic and damping forces are neg- 
lected. 




T^re^hrust force is composed of two components i.e, 



T x = 7*cos(<5 + <D(0)) 



(2-15) 
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T y = Tsin(£ + 0(0)) 



(2-16) 



and a moment applied on the missile is, 

M e = -T cos(<5 + 0(0))v(0) (2-17) 

By applying virtual work principle, the generalized forces of flexible missile is 
found below. 



F> = 



* 



Tcos{6) - T5 sin(0) - TO(O) sin(0) 
T sin(0) + T5 cos(0) + 7*0(0) cos(0) 

-Tv( 0 ) 



(2 - 18) 




ts + rcp(o) 

~Tv( 0) 



(2 - 19) 



where 

T = Force 

<5 = Control Angle being assumed small 
Fy = Large motion generalized force vector 
F u = Small motion generalized force vector 

Appendix C shows the derivation of the generalized forces using the Virtual 
Work Principle. 

The derivation of the equations of motion from the Lagrangian formulation 
yields two sets of coupled equations. One set describes the large motions and the 
another set describes the small motions. These two sets of equations are non- 
linear, coupled, second-order, ordinary differential equations of the form, 



M qq £ + M qn U = F q 



(2 - 20 ) 
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( 2 - 21 ) 



M nq $ + M nn U + G n U + K n U = F n 



where: 

M qq = 3x3 effective inertia matrix for large motions 

M qn = 3x2 coupled inertia matrix of the small motion effect on large 

motions 

F q = 3x1 load vector for the large motions 

M nq = 2x3 coupled inertia matrix of the large motion effect on small 

motions 

M„ n = 2x2 effective inertia matrix for small motions 

G n = 2x2 gyroscopic matrix 

K„ = 2x2 stiffness matrix 

F n = 2x1 load vector for the small motions 

£ = 3x1 vector, generalized coordinates of the large motions 

U = 2x1 vector, generalized coordinates of the small motions 

The detailed development of the equations of motion and definitions of the terms 
are described in Appendix B. 
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III. THE DEVELOPMENT OF A RIGID-BODY CONTROLLER 

In our research, the controller is developed based on the rigid-body assump- 
tion. The purpose of including the rigid-body controller is to perform computer 
simulation and study the dynamic behavior of the flexible missile. The pitch an- 
gle of the missile is controlled by the rigid-body controller which is a single input 
single output (SISO) controller. A general configuration of the flexible missile 
and rigid-body controller are shown in Fig.4. 




figure 4. The flexible missile with rigid-body controller 



The pitch angle is the control variable and the desired states are desired pitch 
angle Q di desired angular velocity 6 d , and desired angular acceleration 6 d . The 
control angle (5) is assumed small. 

The desired error function is defined as 

£ + Kyb 4* Kpt = 6 (3-1) 

where, 
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K p = The position feedback gain. 

K v = The velocity feedback gain. 

e = Q — 6 d = The error in the position. 

i = 8 — Q d = The error in velocity. 

s = 6 — Q d = The error in acceleration. 

The control design begins with the equation of rigid-body (large motion) that 
is obtained by modifying and expanding Eq.(2-20) as 



M n 


M 12 


n i 




fl 


+ 


k 


+ 5 


h 


M 2l 


a/ 22 






fi 




h 2 







Where, 

Vl jAVl.l) M qq (\,2) 

1 11 ['V 2 -') 'V 2 - 2 ). 

M n = [' W w (1,3) 

[M qq (2,3)_ 

M 2i = LM q ft,l) M qq ( 3,2)] 
Mn = lM qq ( 3,3)] 



(3-3) 

(3-4) 

(3-5) 

(3-6) 



h = 



*0 



(3-7) 



?2 = [®] 



(3-8) 



/ = Gravity, centrifugal and coriolis forces. 
f 2 = Gravity, centrifugal and coriolis forces. 
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T cos(0) 
T sin(0) 



h 2 = [0] 



b 



1 - 



-T sin(0) 
T cos(d ) 



b 2 = [0] 



Eq.(3-2) can be rewritten in a tensor form as 



M 11*11 M 12*1 2 ~f\ Jr + b\5 



(3-9) 



(3 - 10) 



(3-11) 



M 2 \*h + M 22 r] 2 =f 2 
From Eq.(3-1 1), 



'/l 




/l + + b\b — Mi 2 >] 2 



Substituting Eq.(3-13) into Eq.(3-12), we Find, 



At] 2 = F — B5 
Where, 

A = M 22 — M 2 \Mi i M12 
F =h ~ M u Mn\fi + *i] 
B = 



(3 - 12) 



(3-13) 



(3 - 14) 



(3 - 15) 

✓ 

(3-16) 

(3-17) 
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Eq.(3-1) can be rewritten as 



(0 - <y + K„(e - <y + k p (6 -e d ) = o 



(3-18) 



By substituting Eq.(3-14) into Eq.(3-18), we find the control angle (5) 



^ = -B-'[A(d d - K v (e - e d ) - K p (e - e d )) - f ] 



(3 - 19) 



The control angle (<5) is a function of the pitch angle (0) , desired pitch angle 
(0 d ), angular velocity (0), desired angular velocity ( 0 d ), desired angular acceler- 
ation ( 0 d ), the position feedback gain Kp , the velocity feedback gain K v , and the 
matrices i.e., A, B, F. Since the thrust force T is included in F and B, the mag- 
nitude of the thrust force will thus affect the control angle <5 . K v and K p are ad- 
justable to obtain the desired response of the pitch angle 0. 
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IV. COMPUTER SIMULATION 



A. SIMULATION OBJECTIVES 

In literature, the application of forces has been limited to gravity, aerodyna- 
mic forces and thrust [Ref. 2]. No damping has been implemented. The thrust 
and aerodynamic forces can be found fairly accurately with standard test and 
design procedures such as static firings to obtain thrust versus time profiles for 
the engine and wind tunnel measurements to determine the aerodynamic forces. 
Gravitational forces can be calculated from the knowledge of the missile's posi- 
tion relative to the earth. The mass including fuel can be estimated from know- 
ledge of the missile's weight before and after burnout (from the measurements), 
and by using a mathematical relationship (often linear) for the decrease in missile 
mass over the engine burn time. 

In this work, the missile's weight is assumed constant and the aerodynamic 
forces are zero. The deformations resulting from the structural flexibility have 
been assumed small and small control angle assumptions are used in the rigid- 
body controller design. 

The primary purpose of this study is to complete the simulation of the flexible 
missile in 2-D motion, where the bending effect is considered to be the only flex- 
ibility source and a rigid-body controller as developed in the last section is in- 
cluded. 

B. SOLUTION TECHNIQUE 

Many numerical integration techniques can be applied in the solution of the 
equations of the motion. The main consideration in the selection of the inte- 
gration technique is the size of the time step necessary to integrate the equations 
of motions with numerical accuracy and stability. 

The type of the equations of motion (2-20, 2-21) in this research permits the 
application of the Sequential Integration Method [Ref. 6], The linear equations 
of small motions are integrated implicitly and the large motion solutions are ob- 
tained using explicit integration method. The implicit methods are effective for 
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linear systems with high frequencies and the explicit methods are effective for 
solving nonlinear systems with low frequencies. The implicit method is especially 
effective for linear systems having a wide range of frequencies of which only the 
lower frequencies are excited. The size of the time step need only be chosen to 
make the solution of the excited modes sufficiently' accurate. Explicit methods 
are effective for large scale systems with low frequencies. Because of the low fre- 
quencies, the size of the time step that we choose for stability need not be small. 
In addition, the explicit method does not need iterative procedures, which are 
time consuming for non-linear systems. 

t 

C. THE COMPUTER SIMULATION CODE 

A high level computer language, i.e., MATLAB, was chosen to simulate the 
missile system. The MATLAB was designed for matrix operations. The simu- 
lation code was developed with modular MATLAB routines. 

The simulation code can be divided into three levels : LEVEL 1 (an overview) 
separates the code into primary portions of INITIALIZATION, PLANT DE- 
SCRIPTION, INTEGRATION, SYSTEM CONTROL and OUTPUT. LEVEL 
2 facilitates several subroutines for LEVEL 1. A listing of MATLAB routines 
required for manipulation in LEVEL 2 is placed in LEVEL 3. 

D. SIMULATION RESULTS 

The computer simulation was performed with variable parameters which de- 
termine missile speed and controller bandwidth. These parameters include force 
T, K v and K p . 

The simulation outputs will be presented in large motions, small motions and 
control angle. The large motions are X Q , Y 0 and 6 and the small motions are v(0) 
and <I>(0). The initial condition for all runs was the pitch angle of 45°. The sim- 
ulation work was divided into two areas. First a simulation was performed for 
the rigid missile using the rigid-body controller. These results were used as a 
baseline for comparison with the results of the flexible missile using the rigid-body 
controller system. Second a simulation was performed for the flexible missile and 
rigid-body controller system. Only the trajectory control will be discussed in the 
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analyze of the dynamics of the flexible missile. The desired trajectory was speci- 
fied as 



0 = 



45° 

45° - 1 12.5r 
0 



when 0.1 < t < 0.1 
when 0.1 ^ t < 0.4 
when t ;> 0.4 



The control angle is assumed to be limited to + /-10 degrees (small angle). 
This restriction puts a saturation line to the control input. 

Fig.5 represents the graphical results of the desired and actual trajectory for 
the rigid missile an'd rigid-body controller. The force (T) and controller band- 
width (co„) were 30000 N and 3 rad/'s respectively. The dark black line and 
dashed line represent the desired and actual trajectory, respectively. The control 
angle (<5) is shown in Fig. 6 and Figures 7-8 present the base positions (A’o.Fo) in 
large motion. 

After presenting results for the rigid missile with rigid-body controller, the 
dynamic behavior of the flexible missile with rigid-body controller will be studied 
next. The force (T) and controller bandwidth (co„) again were 30000 N and 3 
rad/'s, respectively. Fig.9 shows the desired and actual trajectory. The control 
angle is presented in Fig. 10. After 0.1 sec, the missile initially needs a control 
angle which is less than 10 degrees. The effects of small motion can be seen on 
control angle clearly. The base deflection and slope are shown in Figures 11-12. 
There is no small motion effects on flexible missile between 0-0.1 sec because of 
zero control angle and no force components. After 0.1 sec, the small motions are 
excited and have an amplitude of 10 -3 . Figures 13-14 present the base positions 
of the flexible missile in large motion. The large motion behaves like the rigid- 
body motion, which implies that the coupling effect between the large motion and 
small motion is small. The dynamic behavior of the flexible missile is dominant 
by the large motion in this case since the bandwidth of the controller is low and 

S 

small motion is not significantly excited. 
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Figure 5. The desired and actual trajectory for the rigid missile with rigid-body 
controller (T = 30000 N, co n = 3 rad/s) 
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Figure 6. The control angle for the rigid missile with rigid-body controller (T = 
30000 N, co n = 3 rad/s) 
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Figure 7. The large motion position X for the rigid missile with rigid-body controller 
(T = 30000 N, w„ = 3 rad/s) 
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Figure 8. The large motion position Y for the rigid missile with rigid-body controller 
(T = 30000 N, o) n = 3 rad/s) 
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Figure 9. The desired and actual trajectory for the flexible missile with rigid-body 
controller (T = 30000 N, co n = 3 rad/'s) 
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Figure 10. The control angle for the flexible missile with rigid-body controller (T 
= 30000 N, co„ = 3 rad/s) 
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Figure 11. The small motion position v(0) for the flexible missile with rigid-body 
controller (T = 30000 N, co n = 3 rad/s) 
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Figure 13. The large motion position X for the flexible missile with rigid-body con- 
troller (T = 30000 N, co n = 3 rad/s) 
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Figure 14. The large motion position Y for the flexible missile with rigid-body con- 
troller (T = 30000 N, co n = 3 rad/s) 
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Two kinds of tests were performed to further the study of the dynamic be- 
havior of the flexible missile. 

The first test was to study dynamics due to the increase of the control system 
bandwidth (c o„). The desired and actual trajectory of the flexible missile with 
rigid-body controller without saturation on the control angle is presented in 
Fig. 15. The force (T) is still 30000 N while the controller bandwidth (co„) is in- 
creased to 300 rad/'s. As the controller bandwidth is increased, the gap between 
controller bandwidth and system natural frequency is smaller. The system will 
be unstable at controller bandwidths close to the system fundamental frequency 
which can be as high as 270 rad/s in this case. Note that the fundamental fre- 
quency of the simple beam with free-free boundary conditions was calculated as 
co nl =209.2053 rad/s, and the missile's fundamental frequency will be increased 
due to the coupling between large and small motion. The control dynamics in- 
terfere with the structural dynamics. The higher control frequency with the 
bandwidth of 300 rad/s significantly excites the small' motion of the structure. 
The graphs of control angle and the small motions and large motions of the mis- 
sile base (Figures 16-17-18-19-20) show the unstable state. 

The bandwidth of the controller must be set much lower than the funda- 
mental frequency of the missile in order to use the rigid-body control. This implies 
that the pitch-angle response speed, which is determined by the control band- 
width is limited. To keep the response speed, the missile structure must be de- 
signed sufficiently rigid to possess a high fundamental frequency. On the other 
hand, the missile payload will affect the natural frequency of the missile structure 
with the heavier the payload, the lower the fundamental frequency of the missile. 
Therefore, the payload must be limited to achieve high-speed control response. 
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Figure 15. The desired and actual trajectory for the flexible missile with rigid-body 
controller without saturation on the control angle (T = 30000 N, c o n 
. = 300 rad/s) 
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Figure 16. TIfe control angle for the flexible missile with rigid-body controller 
without saturation on the control angle (T = 30000 N, co n = 300 
rad/s) 
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Figure 17. The small motion position v(0) for the flexible missile rigid-body con- 
troller without saturation on the control angle (T = 30000 N, c o n = 
300 rad/s) 
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Figure 18. The small motion position phi(O) for the flexible missile rigid-body con- 
troller without saturation on the control angle (T = 30000 N, co n = 
300 rad/s) 
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Figure 19. The large motion position X for the flexible missile rigid-body controller 
without saturation on the control angle (T = 30000 N, co„ = 300 
rad/s) 
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Figure 20. The large motion position Y for the flexible missile rigid-body controller 
without saturation on the control angle (T = 30000 N, a)„ = 300 
rad/s) 
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Fig. 21 presents the desired and actual trajectory of the flexible missile using 
rigid-body controller with a saturation ( + /-10 Degrees) on the control angle. The 
force (T) and controller bandwidth (c o n ) are 30000 N and 200 rad/s, respectively. 
A better trajectory tracking is obtained, however the actual trajectory is oscillat- 
ing about the desired trajectory which is mainly due to the switching of the con- 
trol angle between the saturation lines. Fig. 22 shows the control angle, and 
Figures 23-24 show the small motion displacement and slope of the base of the 
flexible missile. It was observed that the switching of the control angle between 
+ /-10 degrees causes increased excitation of small motions. Figures 25-26 pre- 
sents the position of missile base in large motion where the large motion is no 
longer dominated by the rigid-body motion. 

The second test was to study the dynamics of the flexible missile due to an 
increase of the missile speed. To increase the speed of the missile, the thrust force 
was increased. Fig. 27 presents the actual and desired trajectory of flexible mis- 
sile. The force was 60000 N and controller bandwidth was 3 rad/s. As seen from 
Fig. 27, the pitch-angle response is determined directly from the bandwidth of the 
controller. The increased missile speed does not change the pattern of the pitch- 
angle response. The control angle of the flexible missile is shown in Fig. 28. The 
control angle is affected by the increased speed. When the speed is increased, the 
control angle gets smaller because a smaller control angle generates a sufficient 
correction to the pitch angle. The increased missile speed has little effect on small 
motions (Figures 29-30). Figures 31-32 show the large motions. 
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Figure 21. The desired and actual trajectory for the flexible missile using rigid-body ' 
controller with saturation on the control angle (T = 30000 N, co„ = 
200 rad/s) 



37 




(saoa)3iDNV noaiNOD 



Figure 22. The control angle for the flexible missile using rigid-body controller with 
saturation on tbe control angle (T =? 30000 N, co„ — 200 rad/'s) 
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Figure 25. The large motion position X for the flexible missile using rigid-body 
controller with saturation on the control angle (T = 30000 N , co n = 
200 rad/s) 
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Figure 26. The large motion position Y for the flexible missile using rigid-body 
controller with saturation on the control angle (T = 30000 N, a ) n 
= 200 rad/s) 
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Figure 27. The desired and actual trajectory for the flexible missile using rigid-body 
controller with increased speed (T = 60000 N, co n = 3 rad/s) 
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Figure 28. The control angle for the flexible missile using rigid-body controller with 
increased speed (T = 60000 N, co n = 3 rad/s) 
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Figure 31. The large motion position X for the flexible missile using rigid-body 
controller with increased speed (T = 60000 N, co n = 3 rad/s) 
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Figure 32. The large motion position Y for the flexible missile using rigid-body 
controller with increased speed (T = 60000 N, co„ = 3 rad/s) 
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V. SUMMARY 



A. CONCLUSIONS 

The development and computer simulation is presented for a bending flexible 
missile with a rigid-body controller system moving in a 2-D coordinate frame. In 
this research, the dynamic model has been developed through the Equivalent 
Rigid Link System utilizing Lagrangian dynamics to obtain a type of system 
equations of motion suited for Sequential Integration Method that integrates 
large motion explicitly and small motion implicitly. The spatial finite element 
discretization of missile structure and the application of truncated natural modal 
responses provide an approximate solution. 

The analysis and simulation were performed to understand the dynamic be- 
havior of a flexible missile using a rigid-body controller. It was found that the 
controller bandwidth must be much lower than the fundamental frequency of the 
missile in order to use the rigid-body controller. The payload will affect the na- 
tural frequency of the missile structure i.e, when the payload is increased, the 
system fundamental frequency will be decreased. The payload must then be lim- 
ited to achieve high-speed response. In order to increase the payload and main- 
tain high-speed control response, a flexible-body controller is needed. 

B. RECOMMENDATIONS 

Areas which remain to be investigated include : 

1. Add the payload and aerodynamic effects to the model. 

2. Design and study the dynamic behavior of a flexible missile with a 
flexible-body controller in 2-D motion. 

3. Build a flexible missile in a laboratory scaie and obtain experimental data. 

4. Design and simulate a control system for flexible missiles in 3-D motion. 
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APPENDIX A. DERIVATION OF THE MODE SHAPE RESPONSE 
MATRIX FOR THE FLEXIBLE MISSILE 

In this study, the natural mode shape functions of a beam are used to repre- 
sent the flexural motion of the flexible missile. Only the first two mode shapes are 
used. The flexible missile is modeled as a continuous Euler-Bernoulli free-free 
beam, neglecting shear deformation and rotary inertia effects. 

The bar (Fig. 33) has system parameters: 

f(x,t) 




Figure 33. The Bar in Flexure 

y(x,t) = Transverse displacement at any point x and time t 

f(x,t) = Transverse force per unit length 

m(x)=The mass per unit length 

EI(x) = Flexural rigidity 

E = Young's modulus of elasticity 

I(x) =The cross-sectional area moment of inertia 

Q(x,t) = Shearing force 

M(x,t) = Bending moment 
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